Drude weight, plasmon dispersion, and a.c. conductivity in doped graphene sheets 
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We demonstrate that the plasmon frequency and Drude weight of the electron liquid in a 
doped graphene sheet are strongly renormalized by electron-electron interactions even in the long- 
wavelength limit. This effect is not captured by the Random Phase Approximation (RPA), com- 
monly used to describe electron fluids and is due to coupling between the center of mass motion 
and the pseudospin degree of freedom of the graphene's massless Dirac fermions. Making use of dia- 
grammatic perturbation theory to first order in the electron-electron interaction, we show that this 
coupling enhances both the plasmon frequency and the Drude weight relative to the RPA value. We 
also show that interactions are responsible for a significant enhancement of the optical conductivity 
at frequencies just above the absorption threshold. Our predictions can be checked by far-infrared 
spectroscopy or inelastic light scattering. 



I. INTRODUCTION 

The first theory of classical collective electron density 
oscillations in ionized gases by Tonks and LangmuiiQ^ in 
the 1920's helped initiate the field of plasma physics. 
The theory of collective electron density oscillations in 
metals, quantum in this case because of higher electron 
densities, was developed by Bohm and PineJ^ in the 
1950's and stands as a similarly pioneering contribution 
to many-electron physics. Bohm and Pines coined the 
term plasmon to describe quantized density oscillations. 
Today plasmonics is a very active subfield of optoelec- 
tronics^ ^, whose aim is to exploit plasmon properties in 
order to compress infrared electromagnetic waves to the 
nanometer scale of modern electronic devices. This wide 
importance of plasmons across different fields of basic and 
applied physics follows from the ubiquity of charged par- 
ticles and from the strength of their long-range Coulomb 
interactions. 

The physical origin of plasmons is very simple. When 
electrons in free space move to screen a charge inho- 
mogeneity, they tend to overshoot the mark. They are 
then pulled back toward the charge disturbance and over- 
shoot again, setting up a weakly damped oscillation. The 
restoring force responsible for the oscillation is the av- 
erage self-consistent field created by all the electrons. 
Because of the long-range nature of the Coulomb in- 
teraction, the frequency of oscillations 0Jpi{q) tends to 
be high and is given in the long wavelength limit by 
oo'^iiq ^ 0) = nq^Vq/m where n is the electron den- 
sity, m is the bare electron mass in vacuum, and Vq is 
the Fourier transform of the Coulomb interaction. This 
simple explicit plasmon energy expression is exact be- 
cause long-wavelength plasmons involve rigid motion of 
the entire plasma, which does not involve the complex 
exchange and correlation effects that dress^ the motion 



of an individual electron. The exact plasmon frequ ency 
expression is correctly captured by the RPA^^^^, but 
also by rigorous argument^ in which the selection of 
a particular center-of-mass position breaks the system's 
Galilean invariance and plasmon excitations play the role 
of Goldstone bosons. In two-dimensional (2D) systems 
Vq = 2Tie^ /q so that uj^\{q ^ 0) = ^J2l\ne^q|m^ where e 
is the magnitude of the electron charge. 

Electrons in a solid, unlike electrons in a plasma or 
electrons with a jellium modef background, experience 
a periodic external potential created by the ions which 
breaks translational invariance and hence also Galilean 
invariance. Solid state effects can lead in general to a 
renormalization of the plasmon frequency, or even to the 
absence of sharp plasmonic excitations. In semiconduc- 
tors and semimetals, however, electron waves can be de- 
scribed at super-atomic length scales using k • p theor}^, 
which is based on an expansion of the crystal's Bloch 
Hamiltonian around band extrema. In the simplest case, 
for example for the conduction band of common cubic 
semiconductors, this leads us back to a Galilean-invariant 
parabolic band continuum model with isolated electron 
energy Edp) = p^/(2mb). The crystal background for 
electron waves appears only via the replacement of the 
bare electron mass by an effective band mass m\^. In this 
type of k p Galilean-invariant interacting electron model, 
valid for many semiconductor and semiconductor hetero- 
j unction systems, the plasmon dispersion is a ccura tely 
given by the random phase approximation ( RPA and 
it is given by the classical formula quoted above with the 
replacements m mb and ^ e^/e, e being the high- 
frequency dielectric constant of the semiconductor mate- 
rial. The absence of electron-electron interaction correc- 
tions to plasmon frequencies at very long wavelengths in 
these systems, has been demonst rated experimentally by 
means of inelastic light scatterin^^^. 
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The situation turns out to be quite different in 
graphene - a monolayer of carbon atoms tightly packed 
in a 2D honeycomb latticd^^^^. When k ■ p theory is 
applied to graphene it leads to a new type of electron 
fluid model, one with separate Dirac-Weyl Hamiltonians 
for electron waves centered in momentum space on one 
of two honeycomb lattice Brillouin-zone corners K and 



\\\\ \ \ \ 



(1) 



Here v is the bare electron velocity, k is the k ■ p 
momentum, a, P are sublattice pseudospin labels, and 
^a(3 = ^ vcctor of Pauli matrices which 

act on the sublattice pseudospin degree-of-freedom. It 
follows that the energy eigenstates for a given p have 
pseudospins oriented either parallel (upper band) or an- 
tiparallel (lower band) to p. Physically, the orientation 
of the pseudospin determines the relative amplitude and 
the relative phase of electron waves on the two distinct 
graphene sublattices. 

Electron-electron interactions in graphene are de- 
scribed by the usual non-relativistic Coulomb Hamilto- 
niaiP 



= Vq Pq p- 



(2) 



where S is the sample area, Vq = 27Te^/{eq) is the 
2D Fourier transform of the Coulomb interaction (e be- 
ing an effective average dielectric constant), and pq = 

Sfc a '^k-q a^k^ci IS the usual density operator. Electron 
carriers with density n can be induced in graphene by 
purely electrostatic means, creating a circular 2D Fermi 
surface in the conduction band with a Fermi radius /cp, 
which is proportional to ^/rl^. The model described by 
1-L = + requires an ultraviolet wavevector cutoff, 
^max, which should be assigned a value corresponding to 
the wavevector range over which describes graphene's 
TT bands. This corresponds to taking k 

max 1/ao where 

ao ~ 1.42 A is the carbon-carbon distance. This model 
is useful when /cmax is much larger than /cp. 

The feature of graphene that is ultimately responsible 
for the large many-body effects on the plasmon dispersion 
and the Drude weight is broken Galilean invariance. The 
lattice reference frame remains present in the continuum 
model through the coupling between momenta and pseu- 
dospins. The oriented pseudospins provide an "ether" 
against which a global boost of the momenta becomes 
detectable. This is explained in detail in the caption of 
Fig.[l] 

Fig. [l] explains why the plasmon frequency in graphene 
is so strongly affected by exchange and correlation. In a 
plasmon mode the region of occupied states (Fermi cir- 
cle) oscillates back and forth in momentum space un- 
der the action of the self-induced electrostatic fleld. In 
graphene, this oscillatory motion is inevitably coupled 
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FIG. 1: Breakdown of Galilean invariance in graphene. Panel 
la) shows the occupied electronic states in the upper band 
of graphene in the ground state. Notice that every state is 
characterized by a value of momentum (the origin of the ar- 
row) and a pseudospin orientation (the direction of the ar- 
row). Panel lb) shows the occupied states after a Galilean 
boost. An observer riding along with the boost would clearly 
see that the orientation of the pseudospins has changed. It 
looks like the pseudospins are subjected to a "pseudomagnetic 
field" that causes them to tilt towards the -\-x direction. The 
appearance of this pseudomagnetic field is the signature of 
broken Galilean invariance. In contrast, in a Galilean invari- 
ant system [Panels Ic) and Id)] the energy eigenstates are 
characterized by momentum only: an observer riding along 
with the boost would not see any change in the character of 
the occupied states. 



with an oscillatory motion of the pseudospins. Since ex- 
change interactions depend on the relative orientation of 
pseudospins they contribute to plasmon kinetic energy 
and renormalize the plasmon frequency even at leading 
order in q. 

In this article we present a many-body theory of this 
subtle pseudospin coupling effect and discuss the main 
implications of our flndings for theories of charge trans- 
port and collective excitations in doped graphene sheets. 
Our manuscript is organized as follows. In Sect. [n| we 
introduce the most important deflnitions and the basic 
linear-response functions that control the plasmon dis- 
persion, the Drude weight, and the a.c. conductivity. In 



Sect. Ill we present the approach we have used to cal- 
culate the Drude weight and the a.c. conductivity, i.e. 
diagrammatic perturbation theory, and the main ana- 
lytical results. Identical results can be obtained using a 
kinetic equation approach which will be detailed in a sep- 
arate publicatioEL' that is focused on a different applica- 
tion. Our main numerical results based on this approach 
are illustrated in Sect. IV In Sect.|V|we emphasize how 
density-density and current-current response functions 
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do not lead to the same results for the Drude weight and 
a.c. conductivity due to the presence of a rigid cutoff in 
momentum space. Finally, in Sect.|VI|we summarize our 



findings and draw our main conclusions. Three Appen- 
dices ( A][C ) highlight some important technical aspects of 
the diagrammatic calculation, while Appendix [P] reports 
on the generalized form of the continuity equation that 
applies in the presence of a rigid momentum cutoff, and 
Appendix [E| reports numerical results for Thomas- Fermi 
screened interactions. 



II. FORMULATION 

The collective (plasmon) modes of the system de- 
scribed by Hamiltonian 1-L can be found by solving the 
following equation^. 



(3) 



where Xpp{q^^) is the so-called proper^^ density-density 
response function. In the q ^ limit of interest here we 
can neglect the distinction between the proper and the 
full causal response function 



XppiQ^^) 
We show below that 



- VqXppiq^^)' 



lim \imjte XppiQ^^) = ^ ^ 



(4) 



(5) 



where P is a as yet unidentified density- and coupling- 
constant-dependent quantity. Note the order of limits in 
Eq. ([5|: the limit ^ is taken in the dynamical sense, 
i.e. vq <^ uj <^ 2eF- Here = '^^f is the Fermi energy 
and 2£f is the threshold for vertical inter-band electron- 
hole excitations. Using Eq. ([5| in Eq. ^ and solving for 
uj we find that, to leading order in 



/ 27re^r 



empi 



-Q , 



(6) 



where we have introduced the "plasmon mass", mpi = 
Tre^n/P. 

In the dynamical limit, the imaginary-part of the a.c. 
conductivity, (j{uj) = ie'^ujliuiq^o Xpp{q^^)/q'^ ^ has the 
form 



V 



(7) 



It then follows from a standard Kramers-Kronig analysis 
that the real-part of the conductivity has a ^-function 
peak at = 0: 5Re g{uj) = V5{uj). Thus the quantity 
V introduced in Eq. ^ is the Drude weight. In the 
presence of disorder the ^-function peak is broadened into 
a Drude peak, but the Drude weight is preserved for weak 
disorder. 



We thus see from Eq. (|6| that the Drude weight com- 
pletely controls the plasmon dispersion at long wave- 
lengths. When electron-electron interactions are ne- 
glected V tends to the RPA Drude weight = 4gF^Q: 
where ctq = e'^/{4h) is the so-called universal^^l^ 
frequency-independent inter-band conductivity of a neu- 
tral graphene sheet. In the same limit mpi hkp/v and 



0) 



gc^e 



(8) 



where g = g^g^ = 4 is a spin-valley degeneracy factor 
and we have introduced the dimensionless fine-structure 
coupling constant o^ee = /{ehv)^ the ratio between the 
Coulomb energy scale e^/cp/e and the kinetic energy scale 
hvkY' (The fine-structure constant can be tuned exper- 
imentally by changing the dielectric environment sur- 
rounding the graphene fiak d^^ l ^^ l) Eq. ^ is the well- 
known RP A^^^^^ ' result for the plasmon dispersion at long 
wavelengths. 

In the following Section we calculate V exactly to 
first order in the fine-structure constant by means 
of diagrammatic perturbation theory, demonstrating in 
the process that its value is substantially enhanced by 
electron-electron interactions. Our results depend on the 
electron density via the ultraviolet cutoff A = /Cmax/^F- 
Thus, the momentum sums that appear in the evalua- 
tion of the diagrams will be restricted in such a way that 
only single-particle states with wave vectors k < /Cmax 
are involved. The value of A varies from ~ 20 for a very 
high-density graphene system with n ~ 10^^ cm~^ to 

100 for a density n ~ 5 x 10^^ cm~^ just large enough 
to screen out the unintended^^ inhomogeneities present 
in samples on substrates. We will see that our results are 
only weakly dependent on A. If this were not true the 
Dirac model for electron-electron interactions in doped 
graphene would not be useful and it would be neces- 
sary to correctly account for interaction effects at energy 
scales beyond those for which the model is valid. The 
cutoff appears only in the well-known electron-electron 
interaction enhancement of the quasiparticle velocity. 



III. DIAGRAMMATIC PERTURBATION 
THEORY 

In Fig. [2] we show the diagrams that contribute to 
the density-density response function Xpp{Qi^) 
first order in the coupling constant; the bare bubble 
diagram, two self-energy-correction diagrams, and one 
vertex-correction diagram. As usual, partial cancella- 
tions between the self-energy and vertex corrections play 
an essential role. In Fig. |2] solid lines are noninteracting 
Green's function^^. 



G(fe,c.) = - E ^ 



(9) 



where U^r is the identity matrix in pseudospin space, 
ak = a k/k, ^fc,^ /J.vk - ep, and r]k^^ r] sgn(-^fc,^) 
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FIG. 2: Feynman diagrams for the proper density- density re- 
sponse function Xpp{Qj^) up to first order in the electron- 
electron interaction. Panel a) The bare bubble diagram. 
Panel b) Vertex correction. Panels c) and d) Self-energy dia- 
grams. 



(with T] = Dashed lines are electron-electron inter- 
actions. 

All wave vectors, /c, and /c', which appear below 
are measured in units of /cp, while frequencies and en- 
ergies are in units of 2eF- Below we set ^ 1. We 
also introduce the density-of-states at the Fermi energy, 

The diagrams in Fig. [2] are first evaluated in the limit 
of small q {i.e. to order q'^) and finite uj. Then, for the 
real part of Xpp we will retain only the terms that scale as 
g^/cj^ for ^ and thus contribute to the Drude weight, 
while for the imaginary part of xpp^ which controls the 
real part of the a.c. conductivity, we calculate the full 
frequency-dependent function. 

Proceeding in this manner we see that the empty bub- 
ble in Fig. [2] {i.e. the noninteracting diagram), de- 
noted by Xpp{Qi reproduces the noninteracting Drude 
weight as expected (see Appendix [A|): 



47rv2g2 ^2 g 



1 / ^Q' 
^{sf)- 

UJ 



(10) 



and 



16 u 
X [e{uj - l)Jv{uJ, A) + {a; ^ -u}] 



(13) 



with 



Jv(w,A) 



V2{uj,l)-Vo{oj,l) 



-V 



1: 



dk 



F{uj,k) 



(14) 



and F{uj, k) = 2fcVi(w, k) + ^[^0(0;, k) + V2{uJ, k)]. Here 
Vm{k,k') are dimensionless Coulomb pseudopotential^, 



V^{k 



exp {—imO) 



(7TF + V + k'^ - 2kk' cos{0) 



(15) 

where gxF = 4aee is the Thomas-Fermi screening wave 
vector (in units of kp). Making use of these formulas it 
is easy to check that the vertex correction to the Drude 
weight [i.e. Eq. (12)] is negative and that the integral in 
the second line of Eq. ( 14) converges in the limit A ^ oo. 



The last two diagrams are "self-energy" corrections, 
denoted by Xpp {Qj'^)j which physically describe the 
modification of the response function due to exchange 
energy corrections to the quasiparticle dispersion. We 
find that (see Appendix [C] for more details) 



1 ^2 



^ dk[Vo{l,k)-V2{l,k)] 



(16) 



and 



and 



X^pJ{q^^) 



.I.j,^s^)^[G{u-l)^{u^-u}] , 

(11) 

where &{x) is the usual Heaviside step function and the 
notation "{cj —a;}" means that we have to add to the 
first term in square brackets an identical term in which 
UJ is interchanged with —uj. 

The next diagram is the so-called "vertex correction" , 
denoted by Xpp\q^^)^ which physically represents the 
dressing of the external driving field by the internally 
generated exchange field. We find that, up to order q'^ 
and for ^ (see Appendix [b| for details). 



where 



16 u 
X [G{uj - 1) Jse(c^, A) + {cj ^ -uj}] 



Jse(c^, A) = -S(c^, A) - a^S(cj, A) 



S(cj,A) 



1 

2 71 



dk kVi{uj^ k) 



(17) 



(18) 



32 

V^2(l,l)] 



Ksf)^[1^o(1,1) 



21^1 (1,1) 
(12) 



In writing Eq. (17) we have excluded a term proportional 



to ^(cj — 1) which is an artifact of perturbation theory as 
we explain below. 
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IV. DRUDE WEIGHT RENORMALIZATION 
AND A.C. CONDUCTIVITY 

We now combine the terms calculated in the previous 
Section. Extensive cancellations occur between vertex 
corrections and self-energy contributions. For example 
the first term on the right hand side of Eq. ( 16 ) cancels 



the vertex contribution (12). The final result for the 



Drude weight to first order in is 



V 



= 1 



(19) 



The real part of the a.c. conductivity is given, to the 
same order, by 



cr(cj) 



= e{uj - 1){1 + aee[Jv{uJ, A) + Jse(c^, A)]} • 

(20) 

Eqs. (19)-(20) are the most important results of this ar- 
ticle. 



A. Long-range interactions 

For unscreened Coulomb interactions [qtf = in 
Eq. ([l5|] ^0(15^) decays as 1/k at large k (see, for ex- 
ample, Ref . l32]) : we thus find that 



V 



where 



= 1 



In (A) 



(21) 



(22) 



is the well-known logarithmic velocity enhancement^^ 
(see also Refs. I32|34|35^) and 



- lim 

4 A^oo 



1 1 



8 47r V3 



dk 
2 



^0(1,^) -^-^2(1,^) 



4C 



In 4 



-0.017 , (23) 



C 0.916 being Catalan's constant. Notice that in 
Eq. ([23| we have taken the limit A ^ 00 since the inte- 
grand decays like l//c^ for /c ^ 00 [/3 reaches its A = 00 
asymptotic value reported in Eq. (23) already at values 
of A as small as ~ 10]. For all dopings of experimental 
relevance A ^ 1: in this regime the doping dependence 
of V is logarithmic and stems from the velocity enhance- 
ment factor (22). The enhancement or suppression of V 



with respect to Vq depends on the relative strength of 
the two terms in Eq. (21), which have opposite sign. In 



the low-density A ^ 1 regime the velocity enhancement 
completely dominates the "/3-term" and V/Vq > 1. The 
enhancement of P, or, equivalently of the plasmon fre- 
quency, can be understood qualitatively by noting that 



electron-electron interactions reduce the pseudospin sus- 
ceptibilit}'^^ and therefore increase the pseudospin stiff- 
ness, i.e. the energy that is required to align pseudospins 
along a given direction. According to the discussion given 
in the Introduction, the larger pseudospin stiffness results 
in higher energy of plasma oscillations. This observation 
is consistent with the fact that in a Landau Fermi liquid 
description^^ the suppression of the pseudospin suscep- 
tibility would be driven by the many-body enhancement 
of the density of states factor v'^ /v^ while interactions 
between the quasiparticles produce the opposite effect. 
Once again we must conclude that the many-body en- 
hancement of the plasmon frequency is intimately con- 
nected to the many-body enhancement of the quasiparti- 
cle velocity v'^ /v. Any physical mechanism that reduces 
v'^ jv without affecting the interaction between quasipar- 
ticles could in principle result in a reduction of the plas- 
mon frequency and Drude weight. 

The Drude weight V and the real part of the a.c. con- 
ductivity G(bj) are plotted as functions of doping n and cj 
respectively in Figs. [3] and [4j We observe that the Drude 
weight (and hence the coefficient of q^l^ in the plasmon 
dispersion relation) is substantially enhanced above the 
noninteracting value. As we showed above in Eqs. ( [21] )- 
(22 ), the enhancement grows slowly (logarithmically) as a 
function of A [it grows linearly with A for short-range in- 
teractions - see Eq. (24) in Sect. 
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. The a.c. conduc- 
tivity is likewise enhanced above the threshold uj — 2£f- 
As we have mentioned above, an unphysical term pro- 
portional to ^{bj — 1) has been omitted from Eq. (20). 
This singular contribution is due to the shift of the inter- 
band absorption threshold from the bare value 25f to the 
dressed value 2£p. In first-order perturbation theory the 
absorption shift appears as a J(cj — 1) contribution to the 
integrand for frequency integrals. 

Fig. [4] shows that for uj > 2£f but uj <C 2A£f the 
a.c. conductivity approaches the high-frequency univer- 
sal value ctq. However, as uj becomes comparable to 2A£:f 
the conductivity decreases and, in fact, has an unphys- 
ical logarithmic divergence at a; = 2A6f. This makes 
perfect sense, since our model is only valid for energies 
that are much smaller that the cutoff energy 2A£:f- As 
long as this condition is met, the calculated spectrum 
is essentially independent of the cutoff. This is a very 
satisfactory feature of the present calculation. 



B. Short-range interactions 

It is instructive to examine how the results presented 
in the main body of this Section change if the electron- 
electron interaction is assumed to be of ultra-short range 
in space, i.e. Vq = v = const and all other moments Vm 
with m > 1 are zero. The calculations can be carried 
out in a completely analytical fashion with the following 
results: 



V 



l + ^(A-l), 



(24) 
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FIG. 3: (Color online) The ratio V/Vq between the interact- 
ing value of the Drude weight P, calculated from Eq. (19), 
and the RPA value Po = 4sFcro is plotted as a function of 
electron density n (in units of 10^^ cm~^) for different values 
of graphene's fine-structure constant ttee- The value aee = 0.9 
is believed to be appropriate for graphene deposited on Si02, 
the other side exposed to air. Note that T> /Vq > 1 and that 
it depends weakly on carrier density. 



where v = ekipv/ {27re'^), and 



3?e (j{uj) 



6(c^-l)[l + aeeJsr(c^,A)] , (25) 



where 



Jsr(cJ, A) 



-1 



In 



(a; + l)(A-cj) 



(cj-l)(A + a;) 



. (26) 



Note that Jsr(co',A) becomes independent of A for A ^ 
oo. It is easy to see that, in this limit, a{uj) approaches 
the universal value for large a;, since Jsr(ct;,oo) goes to 
zero like v/{6uj'^) for ^ oo. If, on the other hand, 
uj is allowed to tend to infinity before A then unphysi- 
cal cutoff-related features appear, such as a logarithmic 
divergence at cj = A. 

Let us make a brief comment on the limit of zero 
doping of our theory. In this limit we predict that the 
strength of the Drude peak vanishes as the absorption 
threshold at 26^ moves toward zero. Then for finite but 
arbitrarily small uj the optical conductivity approaches 
the universal value ctq, yielding a result that is indepen- 
dent of Qfee- 



V. DENSITY RESPONSE VERSUS CURRENT 
RESPONSE 




FIG. 4: (Color online) Panel a) Deviation of the real part of 
the a.c. conductivity cr(uj) from the noninteracting uni- 
versal value (Jo = e^/ (4/i), as a function of frequency uj/{2sf), 
calculated from Eq. ( 20 ) for several values of the fine structure 
constant aee and for n = 1.5 x 10^^ cm~^ (A = 50). Panel 
b) Same as in the main panel but for a fixed value of the 
fine structure constant {aee = 0.9) and two different values of 
doping (corresponding to A = 50 and 100). Notice that the 
dependence of the conductivity on the value of the cutoff A is 
almost invisible at low frequencies and becomes visible only 
at frequencies several times ef- 



In many models of electronic systems, gauge invar i- 
ance and the continuity equation allow us to express the 
density-density response function in terms of the current- 
current response function. In the present model the re- 
lation would take the form 

9 9 

V Q VQ 
XppiQ^^) = -^Xaa{q^^) + ^{[^qJ-q]) , (27) 

where aq is the longitudinal component (parallel to q) 
of the pseudospin-density fluctuation and [aq^p-q] is an 
anomalous commutator^, reminiscent of the commuta- 
tor of Fourier-component-resolved density fluctuations in 
a ID Luttinger liquid^^. Because the current in the 
Dirac model is proportional to the pseudospin density, 
Xaa{Q^ ^) in Eq. (27) is the current-current response func- 
tion. Eq. ( [27| ) works perfectly at the noninteracting level, 
but fails when electron-electron interactions are taken 
into account. Due to the cutoff in momentum space, it 
turns out that the continuity equation. 



idtpq 



Q Jq 



with the current-density operator jq given by 



(28) 



(29) 



is no longer satisfied in the interacting system (see Ap- 
pendix O). 
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Similar conclusions have been reached by 
Mishchenko^^ in the context of calculations of the 
a.c. conductivity of undoped graphene sheets. This topic 
has indeed attracted considerable interest^^ and has 
been at the center of a dispute. Mishchenko^^ , in 
particular, was the first to clarify that the calculation 
of a{uj) based on the density-density response function 
predicts a very small correction over the noninteracting 
conductivity, which is a lso in agreement with the 
experimental findingJ^^^, while methods based on the 
current-current response function^^ or on the kinetic 
equatioii^ predict very different results. We have also 
founcPthat a naive application of Eq. i27l to the inter- 



acting system produces strongly cutoff-dependent results 
for both the Drude weight and the a.c. conductivity of 
the doped system. Moreover, it yields a qualitatively 
different behavior of the Drude weight with respect to 
that found in the present article [see Eq. (19) and Fig.js]. 
The Drude weight calculated from the current-current 
response is suppresse^^ rather than enhanced. On the 
other hand, if Eq. (27) is corrected to take into account 



the modification of the continuity equation due to e-e 
interactions, then the results of the present analysis are 
recovered. The conclusion is that it is generally safer 
in an effective low-energy theory with a rigid cutoff 
to work with the density-density response function 
XppiQ^^) rather than with the current-current response 
function, since the calculation of the former places less 
weight on high-energy intermediate states, which are not 
properly described. This observation is consistent with 
the findings of other author^!. An alternative solution 
is to use a modified Coulomb interaction^^ such that 



the continuity equation (27) remains valid even in the 
interacting system. 



VI. DISCUSSION AND CONCLUSIONS 

Large Fermi velocity enhancements due to exchange 
interactions, like the velocity enhancement that occurs 
in graphene, are common in the theory of solids. Often 
the role of velocity enhancements are fully cancelled in re- 
sponse functions by vertex corrections. For example the 
Fermi velocity of an ordinary two-dimensional electron 
gas diverges when screening is neglected, but the corre- 
sponding reduction in density-density response is absent 
when vertex corrections like those which appear in this 
paper are included. From this point of view the main 
finding of this paper, supported by an explicit first order 
perturbation theory calculation, is that no such cancella- 
tion occurs for graphene's well-known velocity enhance- 
ment. The difference is easy to understand. In the case of 
an ordinary two-dimensional electron gas the velocity en- 
hancement is due to the rapid variation of the unscreened 
exchange self-energy as one goes from occupied states to 
empty states across the Fermi surface. When the density 
changes, the energy at which the velocity peak occurs 
also changes, negating its influence on response. In per- 



turbation theory, this effect is captured by the vertex 
correction. In graphene on the other hand, the velocity 
enhancement occurs over a broad range of energies cen- 
tered on the Dirac point, not the Fermi surface. Density 
response is influenced by enhanced velocities which in- 
crease the energy cost of changing the electron density. In 
ordinary two-dimensional electron gases the logarithmic 
velocity enhancement vanishes in any event once screen- 
ing is accounted for. In graphene, on the other hand 
the enhancement comes from interactions at wavevectors 
much larger than the Fermi energy at which only inter- 
band screening, which does not change the long range 
1/r behavior, is relevant. Neither screening nor vertex 
corrections fully counter the enhanced Drude weight due 
to graphene's Dirac point velocity enhancement. 

We note that a recent experimenlF^ has clearly estab- 
lished the presence of a strong Drude peak which devel- 
ops in graphene as the carrier density is increased. These 
authors conclude that the Drude weight is reduced com- 
pared to the noninteracting electron theory, instead of 
being increased as predicted by this theoretical analysis. 
It will be interesting to see whether or not this exper- 
imental conclusion changes as sample quality improves 
and it becomes possible to more cleanly separate inter- 
band and Drude conductivity contributions over a wider 
regime of carrier density. 

Before concluding we would like to mention that a 
new theoretical paper reporting a study of the effect 
of electron-electron interactions on the conductivity of 
doped graphene steetJ^has appeared recently. (We refer 
the reader to Ref. |47| for studies of band-structure, dis- 
order, phonon, and strain effects.) A quantitative com- 
parison between our findings and theirs is not possible 
since their study takes disorder into account while our 
study is for clean graphene. However, we would like to 
stress that the authors of this paper have treated inter- 
actions only at the vertex level neglecting self-energy ef- 
fects. Even though it is well known that this approxima- 
tion is not "conserving" {e.g. breaks gauge-invariance) 
even in a standard parabolic-band electron liquid, we can 
adopt the same strategy in our calculations and artifi- 
cially "switch-ofF' self-energy contributions to 5Re a{uj) 
[diagrams in panels c) and d) in Fig. |2]. If this proce- 
dure is followed, we find the results shown in Fig. [5] it is 
evident from this point that the neglect of self-energy in- 
sertions is responsible for the reduction of 5Re a{uj) below 
(Jo at large frequencies found in Ref. 46, 

Finally, let us comment on the broader implications of 
our results. Effects similar to those described in this ar- 
ticle occur in graphene bilayers^^ and are also expected 
in other few-layer systems. The lack of Galilean invari- 
ance also affects the cyclotron resonance frequency when 
the 2D sheet of graphene is placed in a perpendicular 
magnetic fielcP^^^ since Kohn's theorenP^, which as- 
serts the absence of many-body effects in cyclotron res- 
onance, is not applicable in this case. (The impact of 
broken Galilean invariance on the collective cyclotron 
motion in graphene has been studied in Ref. [57l in the 
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FIG. 5: (Color online) The real part of the a.c. conductivity 
(in units of ctq) as a function of uj/{2€f) calculated by includ- 
ing only the vertex correction [i.e. by neglecting self-energy 
diagrams in panels c) and d) in Fig. [2] for aee = 0.9 and dif- 
ferent values of doping n (corresponding to A = 100, 200, and 
1000). These data have been calculated using long-range (un- 
screened) interactions, i.e. qtf = in Eq. (15). We remind 
the reader that, in the absence of disorder and within first- 
order perturbation theory, cr(uj) = for cj < 2sf- Notice 
that a^^\uj) converges rapidly in the limit A ^ oo [since, 
as already stated in the main text, the integral in the second 
line of Eq. ( |14[ ) converges in the same limit] and that it drops 
well below the universal value ctq at large cj. 



Appendix A: Calculation of the noninteracting 
contribution to the dynamical density-density 
response function 

The noninteracting response function Xp^p{Qi^)i 
the empty bubble diagram in Fig. [2k), reads 



y(0) 



1 r+oo 1 

{q,co) = gJ2j ^ — Tr[G(fe_,e)G(fc+,e + c.)] 



(27r)2 



E 



IT] 



(Al) 



Here k± = and Tr = ^Tr^r, where Tr^r is the trace 

over pseudospin degrees-of-freedom, and J^pu{k^ k') = 1+ 
jivcos — (l)k')i (j^k being the angle between k and the 
X axis. 

We are interested in the long-wavelength q ^ limit 
of th e de nsity-density response function: we thus expand 
Eq. (Al) to second order in q: 



J),^(/c_,/c+) = (l + /ii^) - jj^iy 



q^ sin^(^fc) 
2/c2 



(A2) 



where we have assumed that q is along the x {i.e. that 
(I) a =_0)- When /i = = — 1, the "Lindhard ratio" in 



Eq. (Al) vanishes. For ji = —u^ Eq. (A2) is already of 
order q"^ and one can just put q = in the Lindhard ratio. 
For fi = u = we need instead to expand the Lindhard 
ratio to second order in q with the result 



high-temperature hydrodynamic regime in which Lan- 
dau levels are not well resolved.) Undoubtedly much 
interesting physics, potentially useful for applications in 
opto-electronics, has still to be learned from the study of 
graphene and othei^ non- Galilean invariant systems. 



IT] 



(7cos((/)fc) S{kF - k) 



vq^ cos^((/)fc)(5 (/cp — k) 



(A3) 
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Here we have used Uk^^^ — n^^^ ^ cos((/)fc)J(/cF — /c)/2, 
and we have dropped irj in the denominator since the 
expansion on the right hand side is always real. The 
first term in Eq. (A3) gives zero contribution to the re- 
spons e fu nction after angula r int egration. Now, replacing 
Eqs. (lA2|) and (lASl) in Eq. (Tni), we find 



Avq^ 



dPk s\v?{(j)k) rik^p - rik^p 



k^ 



2jiivk - 



IT] 

(A4) 



where fi = —fi. Performing the integral over k and in- 
troducing dimensionless variables, one finds Eqs. (10) 



and (11) for the real and imaginary parts of the non- 



interacting response function, respectively. 
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Appendix B: Calculation of the vertex correction 



The vertex correction [i.e. see Fig. [2p)] contribution 
to the density-density response function reads 



1 Z' + OO 7 n-\-00 7 / 

= -s^Jl^^"-"' J ^ij ^Tv[G{k.,e)G{k+,e + co)G{k'+,e' + Lo)G{k'_,e')] 



1 

16 



d2fc f dPk' 



rik' 



rik' 



IT] 



X Tr[(l^ + /io-fc_)(l^ + i/o-fc^)(l^ + i/Vfc^)(l^+/iVfc/_)] . 

I 



(Bl) 



The trace in the previous equation can be expanded in 
powers of q\ 

Tr[(l^ + ll(Tk_){ta + ^CTk^){K + v'(Tk'^){ta + l-l'cTk'J] 

oo 

= 16 ^ </>q; ^, ^0 ^^'^ • 



n=0 



(B2) 



If we choose q = qx/\i\s straightforward to show that 

(i + H(i + /iVO 



ail' (^5 ^ ) 



(B3) 



and 



/ / , ,xSin(0fc)sin(0fc/ - ^fc) 

'^1 = (lJi-y)(lJi -^y ) — ^ — - — 77 ^+Perm , (B4) 

4/c 



where "Perm" is obtained from the first term in Eq. ( B4 ) 



by interchanging primed with non-primed variables. The 



complete expression of W2 is quite cumbersome and will 
not be reported here. The only terms that contribute to 
X^pp{^i^) order are given by 



W2\ 



sin(0fc)sin(0fc/) 



T^^^{k,k') . (B5) 



Collecting all terms up to 0{q^) we can write X^pp^^i^) 
in the following form: 



(B6) 



n=0 



where Xpp denotes terms proportional to Wn 

We find that 



2q^ 



uj ^ v\k- \ — v\kj^ \ ^ IT] uj ^ v\k'_ \ — v\k'j^ \ -\- ir] 

/dPk f dPk' 
(27r)2 J ^^:^^'s-'s'^(^F - /^^)(^(/^F - /^Ocos(^fc)cos(^fc/) [1 + cos(^fc/ - ^fc)] . (B7) 



Introducing the Coulomb pseudopotentials Vm^ performing the integrations over k and k' ^ and using dimensionless 
variable we find immediately Eq. (12). 
For n = 1 we find 



-d[kY - k ) sm(0fc) cos(0fc/) sm(0fc/ - (pk)- 



^ r d'k r dn 

(2^7 (2^ 

p 



)2 k 



UJ + 2/ivk + if] 



Vo{k,kF)-V2{k,kF) 

UJ + 2jiivk + if] 



(B8) 



Since this contribution vanishes for ^ it does not contribute to the renormalization of the Drude weight. Never- 
theless, its imaginary part for finite positive frequency gives precisely the term in the first line of Eq. (14). 
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Finally, for n = 2, we find that 
r.2 r 



2 



— sin((/)fc) sin((/)fc/ 



(1 - nfc,+) (1 - nfc/,+) 



^2 



)+l^2(fc,fcO + 2A^A^'Vi(fc,fcO 
+ iri){u) + 2^'vk' + iry) 

I 



(B9) 



It is possible to show that this expression does not scale 



like 



for 



0: thus it does not contribute to the 



renormalization of the Drude weight. The imaginary part 
of Eq. ( |B9| ) at finite frequency gives precisely the term in 
the second line of Eq. ([l4|. Summing the contributions 

to X^pp{q.i^) coming from the two terms with n = 1 
and n = 2 we find Eq. (13). 



Appendix C: Calculation of the self-energy 
insertions 



We now turn to calculate the two first-order self-energy 
diagrams in Figs.|2]3) and d). The first diagram reads 



+^ de 



Tr 



It is easy to show that for A = the trace becomes 



Tr 



(1^ +/it7fc_)(l^ +i/t7fc^)(l^ +7t7fc^)(l^ + AtTfc^) 

I 



16 J~ niy{k-^^ k—) J- yjy[k-^, k_^) , 



and that the contribution arising from X = —u averages 
to zero after performing the integrations over k and k' . 



Inserting Eq. (C2) in Eq. (CI) we find that 



(Cl) 



(C2) 



^ f cPk f cPk' 



de 



J~ uu{k-\-^ k—)J- yiy(kj^^k'j^) 



2711 (e - + + ^ - + '^Vk+^uY 



^ f d^k f d^k' I AT- /I i'\ 



IT] 



(C3) 



For the diagram in Fig. ^) we find that 



(2^ J (^:^^'^-'^'^'^-'7'^M^(^-'^+)'^7^(^-'^-) 



. (C4) 



r 



We now sum Eqs. (C3) and ( |C4| ) together and separate the inter-band channel is already of order 

inter-band = — /i) from intra-band (y = fi) terms. In 



11 



g^, and thus the other factors can be calculated at g = 0: 



sin^(^fc)cos(^fc/ - ^fc) 



(1 - nfc/,+)(l - nfc,+) 
/c^(co' + 2/ii'/c + if]) 



(C5) 



The real part of this expression vanishes for ^ 0, but 
its imaginary part gives Eq. (17). On the other hand the 



intra-band contribution is entirely real and gives Eq. ( 16 ). 



Appendix D: Generalized continuity equation for an 
interacting system of massless Dirac fermions 

In the presence of a rigid momentum cutoff, the conti- 
nuity equation (27) needs to be modified. To show this we 



start by introducing field operators "^k^a in a restricted 
Hilbert space (RHS): 



(Dl) 



where /cmax is the ultraviolet momentum cutoff and iph.a 
is the regular field operator. All other operators should 
be defined in terms of these new field operators, e.g. the 
density operator reads 



k,a 



(D2) 



Now an interesting observation is that the commutator 
[pq, pq/], which is zero in the regular space, becomes finite 



in the RHS: 

= ^[e(fcn,ax-|fc-g'|)-{g'^q}] 

k,a 

X H-,-,,Jk,a . (D3) 

This is not zero in general. One can show that its expec- 
tation value (. . . )o over the noninteracting ground state 
is zero: ([Pq, Pq'])o = 0- It is also straightforward to show 
that the commutator of the density operator with the ki- 
netic part of the Hamiltonian, remains unchanged 



[pq^y^B] =q jq= Vq&q 



(D4) 



where the Dirac- Weyl Hamiltonian Hd, the current- 
density operator jg, and the longitudinal component of 
the pseudospin-density operator aq are also redefined in 
the RHS. 



As a consequence of Eq. (D3), the commutator of pq 
with the Coulomb interaction He is non-zero in the RHS: 



,p,Q!,/3,7 



^ E E E - |fe - q'\) - {q' ^ q}] H-,-,' ,p^k,a ■ 



q'/O k,a k',13 



(D5) 



Again, one can show that the expectation value of f ^ over spouse function {A = pq^ B = p-q)^ one gets 
the noninteracting ground state vanishes: (tq)o = 0- 
Now applying twice the following identit}^^ 



{{A;B))^ = -{[A,B]) + -{{[A,H];B))^ , (D6) 



where 



{{A 



+ ([rg,/5-g])+l'V(('7g;a_q))^ 

- {(f,;f_,))^ . (D8) 



;B))^ = -i dt ([A(t),5(0)])e-^^^e-^^ (D7) We remind the reader that the linear -response function 
^ Xab{'^) is directly related to the Kubo product by the 

is the usual Kubo product^, to the density-density re- relation Xab{^) = {{^]B))^/S. In the noninteracting 
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limit the first term on the r.h.s. of Eq. (D8) is identically 
zero, the second term is the anomalous commutatoi^^, 
and all terms involving Tq vanish. 

Here , we are mainly interested in the imaginary part 
of Eq. (iDSl): 



2 2 



vq 
1 



(D9) 



The last term on the r.h.s. of Eq. (D9) is at least of 



second order in the electron-electron interaction. Thus, 
up to first order in the Coulomb interaction, we can write 



v'^q^ 

vq 
7a 



(0). 



(DIO) 

where the superscript "(n)" on the response functions 
in the previous equation indicates that they have to be 



evaluated up to the n-th order in the electron-electron 
interaction. 

In order to evaluate xi-r(^'^) Xrv(^'^) start 
from the well known exact-eigestate representation ex- 
pression for a response function at zero temperature^l 



XAB 



{0\A\n){n\B\0) {0\ B \n) {n\ A \0) 



(Dll) 

Before proceeding further let us also introduce the follow- 
ing unitary transformation which diagonalizes the Dirac- 
Weyl Hamiltonian Hb- 



where 



U{k) 



1 /^-i4)kl2 g-i(/)fe/2 



\/2 V 



(D12) 



(D13) 



The full Hamiltonian 1-L after this unitary transformation 
reads 



Qy^O k,k' ,^1,1/,^' ,v' 



(D14) 



r 



where ek,± = ^vk are Dirac-band energies and the ma- 
trix elements /^^/ are defined as 



Now we can write 



(D15) 
. (D16) 



We also introduce 



• (D17) 



xi-r(9,w) 



|a-|n)o(n|f_q|0)o (O|f_q|n)o(n|a-|0); 



= 9 E 



(D18) 



Using the expression of tq given above in Eq. (D5) and Wick's theorem^, and keeping only terms linear in q we 
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finally find that 



r 



_(D191 



In deriving this expression we have assumed that the 
Coulomb interaction is screened {i.e. that Vq is not sin- 
gular for q 0). 



The imaginary part of Eq. (D19) for < < 2vkn 
reads 



^^Xari^^^) = ^ ^ Fp_fc(l - nfc,+)(l - np,+)sin((^p)sin((^p - (/)fc)cos((/)fc)(5(/cmax - /^)(^(^ - 2vp) 



327rv2 



e{uj-2eF)[Vo{uj/{2v),k^,^)-V2{uj/{2v),k^,^)] . 



(D20) 



r 



Using a similar expression for Xrv(^'^) finally obtain 
I 



^rn XariQ^^) - XrliQ^^) = ^^r^^i^ - 2£:F)[Vb(cj/(2v), /c^ax) - 1^2(^/(2v), /cmax)] • 



(D21) 



Taking the limit kmax oo the previous equation sim- 
plifies considerably to 



Appendix E: Numerical results for Thomas-Fermi 
screened interactions 



^ (D22) 



Using this result in Eq. (DIO) we find 



UJ oUJ 



or, in terms of the a.c. conductivity, 



3?e G^^\uj) 



2 2 



(D23) 



-ao . (D24) 



In Figs. [6][8] we present numerical results for P, 
5Re cr(a;), and 3?e (T^^^(cj) obtained by using Thomas- 
Fermi screened interactions, i.e. ^tf 7^ in Eq. (15). 
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